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(57) Abstract: A method for the adaptive adjustment of a PMD compensator in optica] fiber communication systems comprises the 
steps of taking the signal at the compensator output and extracting the components yi(t) and y 2 (t) on the two orthogonal polarizations, 
computing the signal y(t)=[yi(t)j 2 +[y 2 (t)] 2 , sampling the signal y(t) at instants t k =kT with T = symbol interval to obtain samples y(tk), 
computing the mean square error e(k)=y(t k )-u(k) with u(k) equal to the symbol transmitted, and adjusting the parameters of the 
compensator to seek to minimize e(k). A PMD compensator in accordance with the method comprises an adjustment system which 
comprises in turn a photodetecor (17) which takes the components yi(t) and y 2 (t) on the two orthogonal polarizations from the signal' 
at the compensator output, a sampler (19) which samples at instants t*=kT with T = symbol interval, the signal y(t)=[yi(t)] 2 +[y 2 (t)] 2 at 
the output of the photodetector ( 17) to obtain samples y(t0, a circuit (18, 20) for computation of the mean square error e(k)=y(tk)-u(k) 
with u(k) equal to the symbol transmitted, and a regulator (15, 16) which regulates the parameters of the compensator to seek to 
minimize e(k). 
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POLARISATION MODE DISPERSION COMPENSATOR USING A MEAN SQUARE ERROR TECHNIQUE 

The present invention relates to methods of adaptive adjustment of PMD compensators 
in optical fiber communication systems. The present invention also relates to a 
compensator in accordance with said method. 

5 In optical fiber telecommunications equipment the need to compensate the effects of 
polarization mode dispersion (PMD) which occur when an optical signal travels in an 
optical fiber based connection is known. 

It is known that PMD causes distortion and dispersion of optical signals sent over 
10 optical fiber connections making the signals distorted and dispersed. The different time 
delays among the various signal components in the various polarization states acquire 
increasing importance with the increase in transmission speeds. In modern optical fiber 
based transmission systems with ever higher frequencies (10 Gbit/s and more), accurate 
compensation of PMD effects becomes very important and delicate. This compensation 
15 must be dynamic and performed at adequate speed. 

The general purpose of the present invention is to remedy the above mentioned 
shortcomings by making available a method of fast,, accurate adaptive adjustment of a 
PMD compensator and a compensator in accordance with said method. 
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In view of this purpose it was sought to provide in accordance with the present 
invention a method for the adaptive adjustment of a PMD compensator in optical fiber 
communication systems with the compensator comprising a cascade of adjustable 
optical devices over which passes an optical signal to be compensated comprising the 

5 steps of extracting the yl(t) and y2(t) on the two orthogonal polarizations of the signal 
at the compensator output, obtaining the signal y(t)=|yi(t)P+|y2(t)f, sampling the signal 
y(t) at instants t k =kT with T = symbol interval to obtain samples y(t k ), computing the 
mean square error e(k)=y(tk)-u(k) with u(k) equal to the symbol transmitted or with u(k) 
replaced by a decision 0(k) on the transmitted symbol u(k), and producing control 

10 signals for parameters of at least some of said adjustable optical devices to tend toward 
minimization of e(k). 

In accordance with the method it was also sought to realize a PMD compensator in 
optical fiber communication systems applying the method and comprising a cascade of 

15 adjustable optical devices over which passes an optical signal to be compensated and an 
adjustment system comprising a photodetector (17) which takes the components yi(t) 
and y 2 (t) on the two orthogonal polarizations from the signal at the compensator output, 
a sampler 19 which samples at instants tk=kT with T = symbol interval, the signal 
y(t)=|yi( t )P+|y2(0| 2 at the output of the photodetector 17 to obtain samples y(t k ), a circuit 

20 18,20 for computation of the mean square error e(k)=y(tk)-u(k) with u(k) equal to the 
symbol transmitted and a regulator 15,16 which regulates parameters of at least some of 
said optical devices to tend towards minimization of e(k). 
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To clarify the explanation of the innovative principles of the present invention and its 
advantages compared with the prior art there is described below with the aid of the 
annexed drawings a possible embodiment thereof by way of non-limiting example 
applying said principles. In the drawings - 

The components y,(i) and y 2 (t) on the. two orthogonal polarizations are sent 

HG 1 shows a block diagram of a PMD compensator with associated control circuit, 
and 

FIG 2 shows an equivalent model of the PMD compensator. 



With reference to the FIGS FIG 1 shows the structure of a PMD compensator 
designated as a whole by reference number 10. This structure consists of the cascade of 
15 some optical devices which receive the signal from the transmission fiber 11. The first 
optical device is a polarization controller 12 (PC) which allows modification of the 
optical signal polarization at its input. Thus there are three polarization maintaining 
fibers 13 (PMF) separated by two optical rotators 14. 

20 A PMF fiber is a fiber which introduces a predetermined differential group delay 
(DGD) between the components of the optical signal on its two principal states of 
polarization (PSP) termed slow PSP and fast PSP. 
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« 

In the case of the compensator shown in FIG 1 the DGD delays at the frequency of the 
optical carrier introduced by the three PMFs are respectively t c , ate and (1-a) %, with 
0<a<l and with x c and a which are design parameters. 

An optical rotator is a device which can change the polarization of the optical signal 
upon its input by an angle 8i (the figure shows 6, for the first rotator and 6 2 for the 
second) on a maximum circle on the Poincare sphere. 

An optical rotator is implemented in practice by means of a properly controlled PC. 

In FIG 1, xi(t) and x 2 (t) designate the components on the two PSPs of the optical signal 
at the compensator input whereas similarly yi(t) and y 2 (t) are the components of the 
optical signal at the compensator output. 



15 The components y,(t) and y 2 (t) on the two orthogonal polarisations are sent to the input 
of the photo detector (PD) which produces a signal y(t) given by: 

y(t)=|yi(t)| 2 +|y2(t)f 0) 

20 This signal can if necessary be filtered by a post-detection filter. Without loss of 
generality let us assume that this filter is not present. If present this filter produces 
obvious changes in the adaptive adjustment of the compensator parameters. These 
changes are readily imaginable to those skilled in the art and are therefore not further 
described. 
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As shown again in FIG 1 the signal y(t) is then sampled by a sampler 19 at the instants 
tk=kT, where T is the symbol interval. Based on the sample y(t k ) a decision on the 
transmitted bit u(k) is made. Let us designate here by fl(k) this decision made by a 
known decider circuit 18 (DEC). 



The input-output behavior of each of the optical devices is described here by means of 
the so-called Jones transfer matrix H(w) which is a 2x2 matrix characterized by 
frequency-dependent components. Denoting by Wi(co) and W 2 (co) the Fourier 
transforms of the components of the optical signal at the device input the Fourier 
10 transforms Zi(co) and Z 2 (co) of the components of the optical signal at the device output 
are given by: 



(Z X (Q>) 
Z 2 (a» 



W 2 (a» 



(2) 



15 Thus the Jones transfer matrix of the PC is 



-K Kj 



(3) 



20 



where hi e h 2 satisfy the condition |hi| +|h 2 | =1 and are frequency independent. 

Denoting by <|>i and (fe the PC control angles, hi and h 2 are expressed by: 
hi=-cos(<|>2-4 > i) + j sin(4>2-<pi) sin<|»i 
h 2 =-j sin^-fo) cos(|>i 



(4) 



WO 03/050985 PCT/EB02/05661 

6 

Clearly if the PC is controlled using other angles or voltages, different relationships will 
correlate these other parameters with hj and h 2 . The straightforward changes in the 
algorithms for adaptive adjustment of the PMD compensator are discussed below. 

5 Similarly, an optical rotator with rotation angle 9j is characterized by the following 
Jones matrix: 



''cos 6, sin 6,^ 
-sin0, cos0 (> 



(5) 



The Jones transfer matrix of a PMF with DGD n may be expressed as RDR 1 where D is 
10 defined as: 



D = 



0 e- im '' 2 



(6) 



and R is a unitary rotation matrix accounting for the PSPs' orientation. This matrix R 
15 may be taken as the identity matrix I without loss of generality when the PSPs of all the 
PMFs are aligned. 

As shown in FIG 1, to control the PMD compensator a controller 15 is needed to 
produce optical device control signals of the compensator calculated on the basis of the 
20 quantities sent to it by a controller pilot 16 termed controller driver (CD). 

The CD feeds the controller with the quantities needed to update the compensator 
optical device control parameters. As described below, these quantities will be 
extracted by the CD from the signals at the input and/or output of the compensator. 
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The controller will operate following the criterion described below and will use one of 
the two algorithms described below. 

To illustrate the PMD compensator adaptive adjustment algorithms let us assume that 
5 the controller can directly control the parameters fe, fe, 9i and 0 2 which we consolidate 
in a vector 9 defined as: 

If it is not so, in general there will be other parameters to control, for example some 
10 voltages, which will be linked to the previous ones in known relationships. The simple 
changes necessary are discussed below. 

The PMD being a slowly changing phenomenon adjustment of the compensator 
parameters will be performed at a rate lower then the transmitted symbol rate 1/T. Let 
15 us assume that this adjustment is performed at the discrete-time instants tnL=nLT where 
L>1. We designate by: 

9(nL)=(<|>,(nL), k(nL), e «( nL )' 02(nL)) T 
the value of the compensator parameters after the nth update. 

20 

In accordance with the method of the present invention the compensator parameter 
adjusting criterion uses the mean square error (MSE) criterion. 
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Based thereon, the compensator parameters 8 are adjusted to minimize the mean square 
error e(k) defined as: 



e(k)=y(tk)-u(k) 



(7) 



This error is a function of 9 through y(tk). We explicit this dependence by defining 
F(8)=e(k). Therefore, the performance index to be miniinized is the mean value of 
1^(6). The compensator parameters 8 will be updated by the rule: 



10 



9g{F(8)} 



0, [(n + l)L] = ft (nL) - y M^W) = ft (nL) - 2yF (8) 

3 * 8 = 8(«L) 

= MTU L)-2y[y(U)-u(nL)]^^ 



9 = 8(n£.) 



902 



8 = 8(nL) 



9 = 8(nL) 



15 =^(^)-2^)-«(^)]^ 2 



(8) 



ag{F(8)} 



g 1 [(» + l)L] = g,(nL)-r 9£ ^ (9) ^ =^(nL)-2^(8) 

.^i 8 = 8(nL) 1 



8 = 8(nL) 



20 



= 0 2 (nL)-2^(e) 



dE{F(Q)} 



8 = 9(nL) 



30, 



8 = 8(nL) 



^ l(nL) _2^)-« (n L)]^^ 



« * 
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where E{.} denotes "expected" and y>0 is a scale factor which controls the adjustment 
amount. 

In vector notation this means that the vector of the compensator parameters is updated 
5 by adding a new vector with its norm proportional to the norm of the gradient F 2 ^) and 
with opposite direction, i.e. all its components have their sign changed: 

9[(„ + 1)L] = Q(nL) - rVE{F 2 (9)] e = Q{nL) = 6(nL) - 2^(8)VE{F(9)| e = Q(nL) (9) 

10 This way we are sure to move towards a relative minimum of the function F*(9). 

Three variations of the basic updating method defined in (8) can be obtained by using 
only sign information contained in the error e(k) and/or the partial derivative. Hence 
the three possible variations are (considering for example the updating rule related to 

<p 1 [(n+l)L] = q) 1 («L) 2ysign\y{tJ u(nL)] (10) 

( Pl [(n + l)L] = <p 1 (^) 27^) u(nL)]sign . ^ (H) 
20 ( Pl [(« + 1)L] = 9,(^) lysifpiyOm.) u(.nL)hgn ~~ (12) 



We shall now describe two methods of computing the gradient of the function F (9) 
find the required control parameters. 
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First method 



Let us consider the updating rule in vector notation (9). To simplify this rule, in the 
error F(6)=e(k) we substitute the transmitted information symbol u(nL) with the 
5 corresponding decision Q(nL), i.e. we substitute the error e(k) with the estimated error 
e(k) defined as: 

e(k)=y(t k )-Q(k) (13) 

10 In the diagram of FIG 1 this estimated error is obtained at the output of the subtracter 
block 20 and sent to the CD. 

Defining G[9(nL)]=E{e 2 (nL)}, the updating rule (9) becomes: 

15 0 [(* + !) L 1 = 6 (* L ) * Y VG(9)| G = Q(nL) (14) 

The partial derivatives of G(6) for 9=9(t n ) can be computed using the following 5-step 
procedure. 

Step 1: Find the value of G[9(nL)]=G[<j>i(nL), <t> 2 (nL), Gi(nL), 9 2 (nL)] 
20 all'iterazione n. at iteration n. To do this, in the time interval (nLT, nLT+LT/5) 

an estimate of G[9(nL)] is computed by averaging the L/5 values of the 
estimated square error, i.e., 



(15) 
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Step 2: find the partial derivative: 



3G(9) 



9 = 6(nL) 



10 



at iteration n. To do this, parameter fa is set at <|>i(nL)+A while the other 
parameters are left unchanged. The corresponding value of G(9), i.e. 
G[<|>i(nL)+A, <|>2(nL), 8,(nL), 8 2 (nL)] is computed as in step 1 but in the time 
interval (nLT+LT/5, nLT+2LT/5. The estimate of the partial derivative of G(8) 
as a function of <|>i is computed as: 



3G(6) 
Bfa 



9 = 6(nL) 



(nL) + A, fa (nL), 0, (nL), 9 1 (nL)] - G[<t>, (nL), fa (nL), 6, (nL), d 1 (nL)] 

A 

(16) 



15 



Step 3. Find the partial derivative: 



dG(9) 



dfa 



0 = 0(nL) 



20 



at iteration n. To do this the parameter (fe is set at ^(nL)+A while the other 
parameters are left changed. The corresponding value of G(0), i.e. G[(|>i(nL), 
<j>2(nL)+A, G^nL), e 2 (nL)], is computed as in step 1 but in the time interval 
(nLT+2LT/5, nLT+3LT/5). The estimate of the partial derivative of G(6) with 
respect to (fe is computed as: 



3G(6) 
*2 



_ G[fa (nL), fa (nL) + A, 0, (nL), 6 2 (tiL)] - Gfa (nL), fa (nL), g, (nL), 6 1 (nL)) 



0 = Q(nL) 



(17) 
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Step 4: Find the partial derivative: 



dG{9) 



0 = 6(nL) 



10 



at iteration n. To do this, parameter 9i is set at 9i(nL)+ A while the other 
parameters are left unchanged. The corresponding value of G(9), i.e. G[<j>i(nL), 
fcCnL), 8i(nL)+A, erfnL)], is computed as in Step 1 but in the time interval 
(nLT+3LT/5, nLT+4LT/5). The estimate of the partial derivative of G(9) with 
respect to G(9) is computed as: 



3G(9) 



8 = 8(nL) 



Gfa(nLX^(nL),g,(nL) + A.fl a W 

A 

(18) 



15 



Step 5: Find the partial derivative: 



dG(0) 



9 = 9{nL) 



20 



at iteration n. To do this the parameter 9 2 is set at 8 2 (nL)+A while the other 
parameters are left unchanged. The corresponding value of G(9), i.e. G[<t>i(nL), 
MnL), 9i(nL), 9 2 (nL) +A], is computed as in Step 1 but in the time interval 
(nLT+4LT/5, (n+l)LT). The estimate of the partial derivative of G(9) with 
respect to 9 2 is computed as: 



3G(6) 



Gfo(nL),&(iflM(i^),^ 
9 = 9(nL) A (19) 
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The above parameter update is done only after estimation of the gradient has been 
completed. 

Note that in this case it is not necessary that the relationship between the control 
5 parameters of PC and optical rotators and the corresponding Jones matrices be known. 

Indeed, the partial derivatives of the function with respect to the compensator control 
parameters are computed without knowledge of this relationship. Consequentiy if the 
control parameters are different from those assumed as an example and are for example 
10 some voltage or some other angle, we may similarly compute the partial derivative and 
update these different control parameters accordingly. 

Lastly, it is noted that when this algorithm is used the CD is not necessary and the 
controller must receive the estimated error only. 

15 

Second Method 

When an accurate characterization of the PC and each optical rotator is available the 
updating rules may be expressed as a function of the estimated error and the signals on 
20 the two orthogonal polarizations at the compensator input. 

Using the known stochastic gradient algorithm (for example as mentioned in the book 
Digital Communications by J.G.Proakis -McGraw-Hill, New York, 1983) and 
substituting in (8) the error e(nL) with the corresponding estimated error 6(nL) we find: 
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^J(„ + l)L] = ^(«L)-2^(nL)^^ 
<? 2 [(n + 1)L] = ^( n L)-2^(riL)^|^i 



In vector notation the expression (9) becomes: 

9[(« + l)L] = 9(nL) 2y§(nL)Vy(^) 



(20) 



(21) 



10 



Before describing how the gradient of y&L) is to be computed we introduce an 
equivalent model of the PMD compensator. 

It was found that the PMD compensator shown in FIG 1 can be modeled as an 
15 equivalent to a two-dimensional transverse filter using four tapped delay lines (TDL) 
combining the signals on the two principal polarization states (PSP). This equivalent 
model is shown in FIG 2 where: 

c l -cos0 l cosd 2 h i 
c^-sinfljSin&j/ii 

20 c 3 £-sin0 I cos0 2 /£ 

c 4 ^-cosfljsinfl^ 

c s = cos 0j cos 0^ (22) 

Cg^-sin^sin^/ij 

c 7 = sin 0 { cos 9 2 1\ 

Cg^cosfljSintfX 



10 



15 
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For the sake of convenience let c(9)designate the vector whose components are the cj in 
(22). It is noted that the tap coefficients a of the four TDLs are not independent of each 
other. On the contrary, given four of them the others are completely determined by 
(22). In the FIG for the sake of clarity it is designated p=l-a. 

The partial derivatives of y&L) appearing in (21) may be expressed as a function of the 
components on the two PSPs of the signal at the compensator input at some appropriate 
instants. The output sample y(tk) may be written as (where (B) H indicates the 
transposed conjugate of the matrix B): 

y(t k )=c H A(k)c (23) 

where the Hermitian matrix A(k) is given by: 

A(k)=a(k)V(k)+b(k)V(k) (24) 

with vectors a(k) and b(k) defined by: 



20 



a(*) = 



x,{t k -ar c ) 
x x {t k -T c -ar c ) 
*i(f k -ar e ) 

X 2 (h -*c) 

\x 2 it k -t c -az c )j 



Ht k ) = 



' x 2 (t k -2T c ) 

x' 2 (t k -t c ) 
xl{t k -Px c ) 

- x ;(t k -2T e ) 

-x;<t k -T c -0T c ) 

-x\(t t -* e ) 

{ -«t k -(lT e ) ) 



(25) 
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Computing the gradient of y(t n L), the algorithm (9) becomes: 



e[(n+l)L]=9(nL)-4Y e(nL)Re{J H A(nL)c} 



(26) 



5 where 



J = 



'3q 

3ft 
dc 2 

8ft 



9C) 
3& 

3ft 



8c t 
8c 2 

da 



dc 2 

30 2 



(27) 



3c 8 
3ft 



9c 8 



9c g 




10 



is the Jacobean matrix of the transformation c=c(0). 

When the control parameters are different from those taken as examples we will have 
different relationships between these control parameters and coefficients Ci. For 
15 example, if the PC is controlled by means of some voltages, given the relationship 
between these voltages and the coefficients h, and h 2 which appear in (3) by using the 
equations (22) we will always be able to express the coefficients ci as a function of 
these new control parameters. 

20 Consequently in computing the gradient of y(t D L) the only change we have to allow for 
is the expression of the Jacobean matrix J, which has to be changed accordingly as 
readily imaginable to those skilled in the art. 
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Lastly it is noted that when this second method is used the CD must receive the optical 
signals at the input of the compensator and the estimated error d(nL). The CD must 
supply the controller directly with this estimated error and with the signal samples Xi(t) 
and x 2 (t) at the desired instants. 

5 

It is now clear that the predetermined purposes have been achieved by making available 
an effective method for adaptive control of a PMD compensator and a compensator 
applying this method. 

10 Naturally the above description of an embodiment applying the innovative principles of 
the present invention is given by way of non-limiting example of said principles within 
the scope of the exclusive right claimed here. 
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CLAIMS 



1. Method for the adaptive adjustment of a PMD compensator in optical fiber 
communication systems with the compensator comprising a cascade of adjustable 
optical devices over which passes an optical signal to be compensated comprising the 
steps of: 

a. extracting the yi(t) and y 2 (t) on the two orthogonal polarizations of the 
signal at the compensator output, 

b. obtaining the signal y(t)=|yi(t)| 2 +|y 2 (t)| 2 ; 

c. sampling the signal y(t) at instants t k =kT with T = symbol interval to obtain 
samples y(tk), 

d. computing the mean square error e(k)=y(tk)-u(k) with u(k) equal to the 
symbol transmitted or with u(k) replaced by a decision u(k) on the transmitted 
symbol u(k), and 

e. producing control signals for parameters of at least some of said adjustable 
optical devices to tend toward minimization of e(k). 

2. Method in accordance with claim 1 in which said parameters being consolidated 
in a vector 9 and a function F(0)=e(k) being defined the parameters are adjusted to tend 
to minimize the mean value of F 2 ^). 

3. Method in accordance with claim 2 in which the vector 9 being defined as: 

e=(^ 2 ,0 P 0 2 ) T 

where <j>i, <j>2, 9i and 9 2 are said parameters these parameters are updated by the rule: 
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A [(„+l)z,]=^(nL)-yMklwl =^(«L)-27f(e) 
3< *i 9 = 6(nL) 



a£{f(e)} 



9 = 6(nL) 



ag{FQ)} 



8 = 6(/iL) °™ 
^ l(nL) -2r[y(U)-u(nL)]^^ 



dE{F(&)} 



e = 9(nL) 



0 l [( B + l)L] = ^( n Z.)-yM^3 =0,(«L)-2^(e) 

9 = 6(nL) 1 



e=6(/iL) 



a ^ le = 9(nL) 



9 = 6(nL) 



(8) 



where E{.} denotes "expected" and y>0 is a scale factor which controls the adjustment 
amount. 



4. Method in accordance with claim 2 in which the vector 8 of the parameters is 
updated by adding a new vector with the norm proportionate to the norm of the gradient 
of F 2 (9) and with opposite direction, i.e. all its components have their sign changed so 
that the updating rule is: 

9[(n + l)L] = Q(nL) - Y Ve{f 2 (9)} 0 = Q(nL) = W - 2^(9)V£{F(0)| e = Q{nL) 



(9) 



so that movement is towards a relative minimum of the function F 2 ^). 
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5. Method in accordance with claim 2 in which in (8) only the sign of the error 
e(k)=y(t k )-u(k) and/or of the partial derivative is used. 

6. Method in accordance with claim 4 in which in the error F(9)=e(k) the 
transmitted information symbol u(nL) is substituted with the corresponding decision 
fi(nL) so as to substitute the error e(k) with the estimated error e(k) defined as 
e^yCtO-QOc). 

7. Method in accordance with claim 6 in which G[8(nL)]=E{e 2 (nL)} is defined so 
that the updating rule (9) becomes: 

(14) 

e[(» + i)l] = e(«L) - r vg(9)| 9 = q (7iL) 

8. Method in accordance with claim 7 in which the partial derivatives of G(9) for 
0=0(to) are computed by the following five-step procedure: 

- Step 1. Find the value of G[9(nL)]=G[<|>l(nL), <|>2(nL), 91(nL), 92(nL)] at 
iteration n; to do this, in the time interval (nLT, nLT+LT/5) an estimate of 
G[9(nL)] is computed by averaging the L/5 values of the estimated square error, 
i.e., 

L/5 

- Step 2. Find the partial derivative: 



r i XTe(nL+i) 



8G(9) 
d & 9 = Q(nL) 
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at iteration n; to do this, parameter fa is set at <)>i(nL)+A while the other 
parameters are left unchanged; the corresponding value of G(0), i.e. 
G[<|>j(nL)+A, MnL), 9i(nL), 92(nL)] is computed as in step 1 but in the time 
interval (nLT+LT/5, nLT+2LT/5); the estimate of the partial derivative of 
G(8) with respect to <|>i is computed as: 



3G(9) 



8 = 9(nL) 



Gfa (nL) + A, <j> 2 (nL), g, (nL),0 2 (nL)]- G]& (nL), <p 2 (nL), 0, (nL), 6 2 (nL)] 

A 



(16) 



Step 3. Find the partial derivative: 



dG(Q) 



9 = 9(nL) 



at iteration n; to do this the parameter is set at <|)2(nL)+A while the other 
parameters are left unchanged; the corresponding value of G(0), i.e. 
G[<|>i(nL), <|>2(nL)+A, 9i(nL), G 2 (nL)], )], is computed as in step 1 but in the 
time interval (nLT+2LT/5, nLT+3LT/5); the estimate of the partial 
derivative of G(9) with respect to fa is computed as: 



3G(8) 



8 = 9(nL) 



Gfo (nL), & (nL) + A, 6 X {nL), 6 2 (nL)]- Gfo (nL), & (nL), 6 l (nL), d 2 (nL)] 

A 



(17) 



WO 03/050985 



22 



PCT/BB02/05661 



Step 4. Find the partial derivative: 



3G(0) 



30, 



e = 9(nL) 



at iteration n; to do this, parameter 0i is set at 9i(nL)+ A while the other 
parameters are left unchanged; the corresponding value of G(9), i.e. 
G[$i(nL), <(>2(nL), 8i(nL)+A, 9 2 (nL)], )], is computed as in Step 1 but in the 
time interval (nLT+3LT/5, nLT+4LT/5); the estimate of the partial 
derivative of G(9) with respect to G(8) is computed as: 



3G(9) 



80, 



e = 6(nL) 



„ G^, (nL), fa (nL), 0 J (nL) + A, fl 2 (nL)] - G[^, (dL), 0 2 (nL), 6 X (nL), 0 2 (nL)] 

A 

(13) 



Step 5. Find the partial derivative: 

3G(9) 



de i 9 = 9(nL) 

at iteration n; to do this the parameter 9 2 is set at 0 2 (nL)+A while the other 

parameters are left unchanged; the corresponding value of G(8), i.e. • 

G[<|>i(nL), <|>2(nL), 9i(nL), 9 2 (nL) +A], is computed as in Step 1 but in the 
time interval (nLT+4LT/5, (n+l)LT); the estimate of the partial derivative 
of G(9) with respect to 9 2 is computed as: 



8G(9) 



de 2 



6 = 9(nL) 



_ G[(fr, (nL), 0, (nL), 6 X (nL), 6 2 (nL) + A] - gfa (nL), & (nL), g (nL), 0 2 (nL)] 

A 



(.14) 
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9. Method in accordance with claim 3 in which in (8) the error e(nL) is substituted 
with the estimated error 8(nL) and the updating rules become: 

A [(« + l)L] = A W - 2,6(1*0 ^^- } 

^ + i W -^-W^ (20) 

0 2 [(„ + 1)L] = 0 2 (nL) - 2tf*nL) ^^- } 



10. Method in accordance with claim 4 in which the error e(nL) is substituted with 
the estimated error 8(nL) so that the updating rule (9) becomes: 

6[{n + l)L] = 0(nL) - 2^(nL)Vytf, lL ) (21) 

11. Method in accordance with claim 1 in which the PMD compensator is modeled 
as a two-dimensional transverse filter using four tapped delay lines (TDL) combining 
the signals on the two principal polarization states (PSP). 

12. Method in accordance with claims 2 and 11 in which the vector 0 of said 
parameters is updated by adding a new vector with the norm proportionate to the norm 
of the gradient of F 2 (0) and with opposite direction, i.e. all its components have their 
sign changed so that the updating rule is: 

9[(n+l)L]=9(nLHY S(nL)Re{J H A(nL)c} (26) 
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with the Hermitian matrix A(k) given by: 



A(k)=a(k)V(k)+b(k)V(k) (24) 



with the vectors a(k) e b(k) given by: 



a(*) = 



x x {t k -ax c ) 

x 2 (t k ) 
x 2 (t k -ar e ) 
x 2 (t k -r c ) 



b(t k ) = 



f x' 2 (t k -2t e ) N 

xl(t k -T c -pT c ) 

x* 2 (t k -r e ) 

-x\(t k -2x e ) 
-x\{t k -z c -Pt e ) 
-x\if k -r c ) 

- x ;(t k -/k e ) ) 



(25) 



r dc, 


dc x 


3cj 




dA 






de 2 


dc 2 


dc 2 


dc 2 


dc 2 




df 2 


dd l 






dc g 








3& 


30, 


d0 2 

4 



(27) 



with Ci f ....,C8 which are the tap coefficients of the four tapped delay lines and 
xl(t), x2(t) are the components of the two principal polarization states at the 
compensator input. 
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13. Method in accordance with claim 1 in which said optical devices comprise a 
polarization controller with control angles fc, and two optical rotators with rotation 
angles 0i and 9 2 and said parameters comprise said control angles <J)i, <fe and said 
rotation angles 0i, 62 or functions thereof. 

14. Method in accordance with claim 13 in which between the controller and an 
optical rotator and between optical rotators there are fibers which introduce a 
predetermined differential unit delay maintaining the polarization. 

15. Compensator for PMD in optical fiber communication systems applying the 
method in accordance with any one of the above claims and comprising a cascade of 
adjustable optical devices over which passes an optical signal to be compensated and an 
adjustment system comprising a photodetector (17) which takes the components yi(t) 
and y 2 (t) on the two orthogonal polarizations from the signal at the compensator output, 
a sampler (19) which samples at instants t k =kT with T = symbol interval, the signal 
y(t)=|yi(t)| 2 +|y2(0| 2 at the output of the photodetector (17) to obtain samples y(tk), a 
circuit (18,20) for computation of the mean square error e(k)=y(tk)-u(k) with u(k) equal 
to the symbol transmitted and a regulator (15,16) which regulates parameters of at least 
some of said optical devices to tend towards minimization of e(k). 

16. Compensator in accordance with claim 15 characterized in that said optical 
devices comprise a polarization controller with control angles <j>i, <k and two optical 
rotators with rotation angles 61 and 9 2 and in which said parameters which are adjusted 
consist of said control angles $u fa and said rotation angles 0i and 8 2 . 
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17. Compensator in accordance with claim 16 characterized in that between the 
controller and optical rotator and between optical rotators there are fibers which 
introduce a predetermined differential unit delay maintaining the polarization. 
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